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Abstract 

We study the spin-| Hubbard-Kondo lattice model by means of the Composite Operator Method, after applying 
a Holstein-Primakov transformation. The spin and particle dynamics in the ferromagnetic state are calculated by 
taking into account strong on-site correlations between electrons and antiferromagnetic exchange among | spins, 
together with usual Hund coupling between electrons and spins. 
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The revival in the study of manganites has led to 
experimental re-examination of their different proper- 
ties. One of the puzzling features is the strong deviation 
of the spin-wave dispersion from the typical Heisen- 
berg behavior. In particular, it has been observed an 
unexpected softening at the zone boundary [1]. These 
observations are very important as they indicate that 
some aspects of spin dynamics in manganites have not 
been entirely understood yet [2,3]. According to this, 
we have decided to investigate the spin dynamics of 
the ferromagnetic state of the spin-| Ifubbard-Kondo 
lattice model by means of the Composite Operator 
Method (COM) [4]. The Hamiltonian under analysis 
reads as 

H^Y. i-^dta^i - fiSii) (i) c (j) + Uj2n^ (i) (i) 

ij i 

- Jh ^ S (j) ■ S (i) + dJAF ^ S (i) • (i) (1) 
i i 

where i is a vector of the d-dimensional lattice and i — 
(i, t), /i is the chemical potential, c' (i) = 
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is the electronic creation operator in spinorial nota- 
tion, t is the hopping amplitude, Qij is the nearest- 
neighbor projector, U is the on-site Coulomb interac- 
tion, rio- (i) = cl (i) Co- (i), Jh is the Hund coupling, 
s (i) = ic^ (i) ere (i), a are the Pauli matrices, S («) is 
the core |-spin, Jaf is the antiferromagnetic coupling. 
We have used the notation (j>°' (i) = X^j Qij^ (i, t)- 

In order to avoid the difficulties related to the high 
value of the core spin, we have used the Holstein- 
Primakoff transformation: S+ (i) — \/2S — Ua {i)a (i), 
(i) = flt (i) ^2S-na (i), (i) = S-na (i), where 
S = I , a (i) is a spinless bosonic destruction operator 
and Ua (i) — a^{i)a{i). Then, we have decided to 
approximate the non-linear term ^2S — Ua (i) to the 
first order in 2S-^n where Sua (i) = ria {i) ~ {ria (i)} 
and Ua — {ua (i)). It is worth noticing that this ap- 
proximation preserves all properties related to the 
angular momentum algebra of the core spin. The 
transformed Hamiltonian reads as 

H = Y1 (-2rftoij - Mf^ij) c+ {i) c (j) + UY,n^ (i) (i) 

ij ■ 
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- 2dJAF {S - ria) (^) 

i 

+ dJAFvl^^ot(i)o"(i) (2) 



where A = \/2S — Ua- 

Within the framework of the COM, we have chosen 
two opcratorial basis to study the spin and particle 

dynamics 



B(i) = 



a{i) 
s+ (i) 



( (i) \ 



(3) 



where ^ (i) = [1 — n c (i) and ri{i) = n (i) c (i). 

Then, we have hnearized the equations of motion by 
projecting the currents onto the basis 



where 



3 

eB,F (i, j) = mB,F (i, 1) (IJ) 



rriB (ij) 

mp (iJ) 
Ib (i, j) 



i|v>(i,t)^Uj,^)}) 
[b {i,t)BHi,t)] 



(4) 
(5) 

(6) 

(7) 

(8) 
(9) 



/F(iJ) = ({^(i,i)VUj,t)}) (10) 

This procedure assures that the neglected component 
of the current is orthogonal to the chosen basis. Ac- 
cording to this, we have obtained the corresponding 
retarded Green's functions in the pole approximation 



' B,F 



(k) 



(k) + i^ 



(11) 



where the energies -E^*'^ (k) are the eigenvalues of the 
energy matrices £s,f (k) = J^[eB,F (i, j)] and the spec- 
tral densities Ug^p (k) can be computed in terms of the 

normalization matrices Ib,f and of the eigenvectors of 
the energy matrices[4]. is the Fourier transform. 

The parameters appearing in the expressions of 
mB,F and Ib,f have the following definitions 



n = {n (i)> = 2 - {Cfh + Cf22 + Cf33 + Cfh) (12) 



m = (s^ (i)) = - (Cf44 - Cf22) 

"a = CbII — 1 

Pi — ( C'fII + ^F22 + ^F33 + ^"^744 ) 

1 (a(i) s_ (j)> _ 1 Cbi2 



P2 



P5 ■ 



A m -1 m 

1 (fa„(0.sa»)) 
m 

1 ((Sria (i) a (i) s_ (i)) 
1 (i5n„ (i) n (i)) 



F33 



^0—2 ~ '-'F22 + 

xl = (n° («) (i)) 



^"44 / 



F22 '-'P33 + ^F44j 

s° (i) s(i)^^ 



(13) 

(14) 
(15) 

(16) 
(17) 
(18) 
(19) 
(20) 
(21) 

(22) 
(23) 



where Cfi,.^ = (^/3 (i) (i)), C^f^-, = (i) i>\ (i)> 
and CBfi-i ~ (5/3 (i) (i)) • According to the prescrip- 
tions of the COM [4], the parameters that cannot be 
computed by their definitions (ps, p4, ps, p and Xz) 
would be fixed by the following relations 



Cii = C33 C12 = 
C34 = C44 = Cbii 



(24) 
(25) 



coming from the algebra and by one more relation 
coming from the request that the hydrodynamic limit 
should be satisfied (i.e., the existence of a sound mode). 
All these equations (definitions and constraints) form 
a coupled system that should be computed self- 
consistently. The results of these calculations will be 
presented elsewhere. 

In conclusion, we have reported the solution for 
the Hubbard-Kondo model in presence of antiferro- 
magnetic coupling between the core spin with the 
framework of the Composite Operator Method in the 
pole approximation. The model has been first mapped 
through the Holstein-Primakov transformation that 
has been then approximated to the first order in the 
number fluctuation operator. The spin dynamics has 
been fully determined and will be analyzed numeri- 
cally. 



References 

[1] H. Y. Hwang, et al., Phys. Rev. Lett. 80 (1998) 1316. 

[2] G. KhaliuUin, P. Kilian, Phys. Rev. B 61 (2000) 3494. 

[3] F. Mancini, N. B. Perkins, N. M. Plakida, Phys. Lett. A 
284 (2001) 286. 



2 



[4] F. Mancini, A. Avella, Adv. Phys. 53 (2004) 537. 



3 



